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More than half of all stars in the local Universe are found in massive 
spheroidal galaxies 1 , which are characterized by old stellar populations 2 ' 3 
with little or no current star formation. In present models, such galaxies 
appear rather late as the culmination of a hierarchical merging process, in 
which larger galaxies are assembled through mergers of smaller precursor 
galaxies. But observations have not yet established how, or even when, the 
massive spheroidals formed 2 ' 3 , nor if their seemingly sudden appearance 
when the Universe was about half its present age (at redshift z rs 1) results 
from a real evolutionary effect (such as a peak of mergers) or from the 
observational difficulty of identifying them at earlier epochs. Here we 
report the spectroscopic and morphological identification of four old, fully 
assembled, massive (> 10 11 solar masses) spheroidal galaxies at 1.6 < z < 1.9, 
the most distant such objects currently known. The existence of such 
systems when the Universe was only one-quarter of its present age, shows 
that the build-up of massive early-type galaxies was much faster in the 
early Universe than has been expected from theoretical simulations 4 . 

In the ACDM scenario 5 , galaxies are thought to build-up their present-day mass 
through a continuous assembly driven by the hierarchical merging of dark matter 
halos, with the most massive galaxies being the last to form. However, the formation 
and evolution of massive spheroidal early-type galaxies is still an open question. 

Recent results indicate that early-type galaxies are found up to z ~ 1 with a num- 
ber density comparable to that of local luminous E/SO galaxies 6 ' 7 , suggesting a slow 
evolution of their stellar mass density from z ~ 1 to the present epoch. The critical 
question is whether these galaxies do exist in substantial number 8 ' 9 at earlier epochs, 
or if they were assembled later 10,11 as favored by most renditions of the hierarchical 
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galaxy formation scenario 4 . The problem is complicated also by the difficulty of iden- 
tifying such galaxies due to their faintness and, for z > 1.3, the lack of strong spectral 
features in optical spectra, placing them among the most difficult targets even for the 
largest optical telescopes. For example, while star-forming galaxies are now routinely 
found up to z ~ 6.6 12 , the most distant spectroscopically confirmed old spheroid is 
still a radio-selected object at z — 1.552 discovered almost a decade ago 13 ' 14 . 

One way of addressing the critical question of massive galaxy formation is to search 
for the farthest and oldest galaxies with masses comparable to the most massive 
galaxies in the present-day universe (10 11-12 M Q ), and to use them as the "fossil" 
tracers of the most remote events of galaxy formation. As the rest-frame optical - 
near-infrared luminosity traces the galaxy mass 15 , the K s -band (A ~ 2.2 /im in the 
observer frame) allows a fair selection of galaxies according to their mass up to z ~ 2. 

Following this approach, we recently conducted the K20 survey 16 with the Very 
Large Telescope (VLT) of the European Southern Observatory (ESO). Deep optical 
spectroscopy was obtained for a sample of 546 objects with K s < 20 (Vega photomet- 
ric scale) and extracted from an area of 52 arcmin 2 , including 32 arcmin 2 within the 
GOODS-South field 17 (hereafter the GOODS/K20 field). The spectroscopic redshift 
(zspec) completeness of the K20 survey is 92%, while the available multi-band pho- 
tometry (BVRIzJHK s ) allowed us to derive the spectral energy distribution (SED) 
and photometric redshift (z phot ) of each galaxy. The K20 survey spectroscopy was 
complemented with the ESO/GOODS public spectroscopy (Supplementary Table 1). 

The available spectra within the GOODS/K20 field were then used to search for 
old, massive galaxies at z > 1.5. We spectroscopically identified four galaxies with 
18 < K s < 19 and 1.6 < z spec < 1.9 which have rest-frame mid-UV spectra with 
shapes and continuum breaks compatible with being dominated by old stars and R — 
K s > 6 (the colour expected at z > 1.5 for old passively evolving galaxies due to the 
combination of old stellar populations and k-correction effects 9 ). The Supplementary 
Table 1 lists the main galaxy information. The spectrum of each individual object 
allows a fairly precise determination of the redshift based on absorption features and 
on the overall spectral shape (Fig. 1). 

The co-added average spectrum of the four galaxies (Fig. 2-3) shows a near-UV 
continuum shape, breaks and absorption lines that are intermediate between those 
of a F2 V and a F5 V star 18 , and typical of about 1-2 Gyr old synthetic stellar 
populations 19,20 . It is also very similar to the average spectrum of z ~ 1 old Extremely 
Red Objects 7 (EROs), and slightly bluer than that of the z ~ 0.5 SDSS red luminous 
galaxies 21 and of the z = 1.55 old galaxy LBDS 53w091 13 . However, it is different in 
shape and slope from the average spectrum of z ~ 1 dusty star-forming EROs 7 . 

The multi-band photometric SED of each galaxy was successfully fitted without 
the need for dust extinction, and using a library of simple stellar population (SSP) 
models 19 with a wide range of ages, Z = Z Q and Salpeter IMF. This procedure yielded 
best-fitting ages of 1.0-1.7 Gyr, the mass-to-light ratios and hence the stellar mass of 
each galaxy, which results in the range of 1-3 xlO 11 hj 2 M . H = 70 km s" 1 Mpc -1 
(with h-jo = Ho/ 70), Q m = 0.3 and Q\ = 0.7 are adopted. 

In addition to spectroscopy, the nature of these galaxies was investigated with the 
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fundamental complement of Hubble Space Telescope+ ACS (Advanced Camera for 
Surveys) imaging from the GOODS public Treasury Program 17 . The analysis of the 
ACS high-resolution images reveals that the surface brightness distribution of these 
galaxies is typical of elliptical/early- type galaxies (Fig. 4). 

Besides pushing to z ~ 1.9 the identification of the highest redshift elliptical 
galaxy, these objects are very relevant to understand the evolution of galaxies in 
general for three main reasons: their old age, their high mass, and their substantial 
number density. 

Indeed, an average age of about 1-2 Gyr (Z = Z & ) at <z>~ 1.7 implies that the 
onset of the star formation occurred not later than at z ~ 2.5 — 3.4 (z ~ 2 — 2.5 for Z = 
2.5Z Q ). These are strict lower limits because they follow from assuming instantaneous 
bursts, whereas a more realistic, prolonged star formation activity would push the bulk 
of their star formation to an earlier cosmic epoch. As an illustrative example, the 
photometric SED of ID 646 (z = 1.903) can be reproduced (without dust) with either 
a ~1 Gyr old instantaneous burst occurred at z ~ 2.7, or with a ~2 Gyr old stellar 
population with a star formation rate declining with expi—tjr) (r = 0.3 Gyr). In the 
latter case, the star formation onset would be pushed to z ~ 4 and half of the stars 
would be formed at z ~ 3.6. In addition, with stellar masses M* > lO n /if o 2 M , these 
systems would rank among the most massive galaxies in the present-day universe, 
suggesting that they were fully assembled already at this early epoch. 

Finally, their number density is considerably high. Within the comoving volume 
relative to 32 arcmin and 1.5 < z < 1.9 (40,000 Mpc 3 ), the comoving density of 
such galaxies is about 10 -4 h^ Mpc -3 , corresponding to a stellar mass density of about 
2 x 10 7 hjo M Mpc -3 , i.e. about 10% of the local [z = 0) value 22 for masses greater 
than 10 11 M . This mass density is comparable to that of star-forming M* > lO n M 
galaxies at z ~ 2 23 , suggesting that while the most massive galaxies in the local 
universe are now old objects with no or weak star formation, by z ~ 2 passive and 
active star-forming massive galaxies coexist in nearly equal number. 

Although more successful than previous models, the most recent realizations of 
semi-analytic hierarchical merging simulations still severely underpredict the density 
of such old galaxies: just one old galaxy with K s < 20, R — K s > 6, and z > 1.5 is 
present in the mock catalog 4 for the whole five times wider GOODS/CDFS area. 

As expected for early- type galaxies 9,24 , the three galaxies at z ~ 1.61 may trace 
the underlying large scale structure. In this case, our estimated number density may 
be somewhat biased toward a high value. On the other hand, the number of such 
galaxies in our sample is likely to be a lower limit due to the spectroscopic redshift 
incompleteness. There are indeed up to three more candidate old galaxies in the 
GOODS/K20 sample with 18.5 < K s < 19.5, 1.5 < z phot < 2.0, 5.6 < R - K s < 6.8 
and compact HST morphology. Thus, in the GOODS/K20 sample the fraction of old 
galaxies among the whole z > 1.5 galaxy population is 15±8% (spectroscopic redshifts 
only), or up to 25=1=11% if also all the 3 additional candidates are counted. 

It is generally thought that the so-called "redshift desert" (i.e. around 1.4 < z < 
2.5) represents the cosmic epoch when most star formation activity and galaxy mass 
assembly took place 25 . Our results show that, in addition to actively star forming 
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galaxies 26 , also a substantial number of "fossil" systems already populate this redshift 
range, and hence remain undetected in surveys biased towards star-forming systems. 
The luminous star-forming galaxies found at z > 2 in sub-mm 27 and near-infrared 23 ' 28 
surveys may represent the progenitors of these old and massive systems. 
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Figure 1 

The individual and average spectra of the detected galaxies. From bottom to top: 
the individual spectra smoothed to a 16 A boxcar (26 A for ID 237) and the average 
spectrum of the four old galaxies (z average = 1.68). The red line is the spectrum of 
the old galaxy LBDS 53w091 (z = 1.55) used to search for spectra with a similar 
continuum shape. Weak features in individual spectra (e.g. MgIIA2800 and the 
2640 A continuum break, B2640) become clearly visible in the average spectrum. 
The object ID 235 has also a weak [OII]A3727 emission (not shown here). The spectra 
were obtained with ESO VLT+FORS2, grisms 2001 (R(l") ~400) (ID 237) and 3001 
(R(l") ~600) (IDs 235,270,646), 1.0" wide slit and < 1" seeing conditions. The 
integrations times were 3 hours for ID 237, 7.8 hours for IDs 235 and 270, and 15.8 
hours for ID 646. For ID 646, the ESO/GOODS public spectrum was co-added to 
our K20 spectrum (see Supplementary Tab. 1). "Dithering" of the targets along the 
slits was applied to remove efficiently the CCD fringing pattern and the strong OH 
sky lines in the red. The data reduction was done with the IRAF software package 
(see 16 ). The spectrophotometric calibration of all spectra was achieved and verified by 
observing several standard stars. The average spectrum, corresponding to 34.4 hours 
integration time, was obtained by co-adding the individual spectra convolved to the 
same resolution, scaled to the same arbitrary flux (i.e. with each spectrum having the 
same weight in the co- addition), and assigning wavelength-dependent weights which 
take into account the noise in the individual spectra due to the OH emission sky lines. 

Figure 2 

The detailed average spectrum of the detected galaxies. A zoom on the average 
spectrum (blue) compared with the synthetic spectrum 19 of a 1.1 Gyr old simple 
stellar population (SSP) with solar metallicity (Z = Z Q ) and Salpeter IMF (red). 
The observed average spectrum was compared to a library of synthetic SSP template 
spectra 19 ' 20 with a range of ages of 0.1-3.0 Gyr with a step of 0.1 Gyr, and with assumed 
metallicities Z=0Ax, l.Ox, and 2.5xZ & . The best fit age for each set of synthetic 
templates was derived through a \ 2 minimization over the rest-frame wavelength 
range 2300-3400 A. The rms as a function of wavelength used in the \ 2 procedure 
was estimated from the average spectrum computing a running mean rms with a step 
of 1 A and a box size of 20 A, corresponding to about three times the resolution of the 
observed average spectrum. The median signal-to-noise ratio is ~20 per resolution 
element in the 2300-3400 A range. The wavelength ranges including the strongest 
real features (i.e. absorptions and continuum breaks) were not used in the estimate 
of the rms. The resulting reduced x 2 is of the order of unity for the best fit models. 
In the case of solar metallicity, the ranges of ages acceptable at 95% confidence level 
are I.OIqi Gyr and 1.4l°,4 Gyr for SSP models of 19 and 20 respectively (see also Fig. 
3, top panel). Ages ~ 50% younger or older are also acceptable for Z = 2.5Z Q or 
Z = 0AZ Q respectively. The 2640 A and 2900 A continuum break 13 amplitudes 
measured on the average spectrum are B2640=1.8±0.1 and B2900=1.2±0.1. These 
values are consistent with the ones expected in SSP models 19-20 for ages around 1-1.5 
Gyr and solar metallicity. For instance, the SSP model spectrum shown here has 
B2640=1.84 and B2900=1.27. 
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Figure 3 

The comparison between the average spectrum and a set of spectral templates. The 
average spectrum (blue) compared to a set of template spectra. From bottom: F2 V 
(green) and F5 V (red) stellar spectra 18 with Z = Z Q , the composite spectrum (red) 
of 726 luminous red galaxies at 0.47 < z < 0.55 selected from the SDSS 21 (available 
only for A > 2600 A), the average spectra of z ~ 1 old (red) and dusty star- forming 
(green) EROs 7 , SSP synthetic spectra 19 (Z = Z Q , Salpeter IMF) with ages of 0.5 Gyr 
(magenta), 1.1 Gyr (green) and 3.0 Gyr (red). 

Figure 4 

The morphological properties of the detected galaxies. Images of the four galax- 
ies taken with the Hubble Space Telescope +ACS through the F850LP filter (from 
GOODS data 17 ) which samples the rest-frame ~3000-3500 A for 1.6 < z < 2. The im- 
ages are in logarithmic grey-scale and their size is 2" x 2", corresponding to ~ 17 x 17 
kpc for the average redshift z — 1.7 and the adopted cosmology. At a visual inspection, 
the galaxies show rather compact morphologies with most of the flux coming from the 
central regions. A fit of their surface brightness profiles was performed with a "Sersic 
law" (oc r l l n ) convolved with the average point spread function extracted from the 
stars in the ACS field and using the GASPHOT 29 and GALFIT 30 software packages. 
Objects ID 237 and ID 646 have profiles with acceptable values of n in the range of 
4 < n < 6, i.e., typical of elliptical galaxies, object ID 270 is better reproduced by 
a flatter profile (1 < n < 2), whereas a more ambiguos result is found for the object 
showing some evidence of irregularities in the morphology (ID 235, 1 < n < 3). These 
latter objects may be bulge-dominated spirals but no bulge/disk decomposition was 
attempted. Ground-based near-infrared images taken under 0.5" seeing conditions 
with the ESO VLT+ISAAC through the K s filter (rest-frame ~6000-8000 A) show 
very compact morphologies, but no surface brightness fitting was done. 
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SUPPLEMENTARY TABLE 1 
The detected galaxies 



IAU K20 R.A. (J2000) Dec (J2000) K s R-K s z Spectrum 

ID ID h m s o I II 

J033210. 79-274627.8 235 03 32 10.776 -27 46 27.73 17.98±0.04 6.47±0.10 1.610 K20 

J033210. 52-274628.9 237 03 32 10.507 -27 46 28.84 19.05±0.05 6.83±0.28 1.615 K20 

J033212. 53-274629.2 270 03 32 12.525 -27 46 29.16 18.74±0.05 5.99±0.10 1.605 K20 

J033233.85-274600.2 646 03 32 33.847 -27 46 00.24 19.07±0.07 5.99±0.10 1.903 K20+GOODS 



Supplementary Table 1 

IAU ID: official identification number in the GOODS-South catalog (z-band) 



(http:/ /www.stsci.edu/science/goods/catalogs). 



K20 ID: identification number in the K20 survey catalog ( |http : / / www. arcetri . astro . it / ^k2 Wf\ . 

R.A., Dec: Right Ascension and Declination at equinox J2000 based on the pub- 
lic ESO/GOODS AVband VLT+ISAAC image. 

K s : K20 survey total magnitude in the AVband (Vega scale). 

R — K s color (Vega scale) in 2" diameter aperture. 

z: spectroscopic redstart . 

Spectrum: K20: K20 survey, GOODS: public ESO/GOODS VLT+FORS2 spec- 



troscopy (Vanzella et al., in preparation; |http:/ /www.eso.org/science/goods ). 
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Figure 2: 
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